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ABSTRACT 


This  thesis  consists  of  two  parts.  First,  an 
experimental  investigation  of  a  new  device  called  the 
whistler  nozzle  was  conducted.  Experiments  were 
conducted  in  the  areas  of  nozzle  efficiency,  mass 
entrainment,  and  flow  visualization.  Flow 
visualization  shoved  the  presence  of  a  Coanda  type  jet 
wall  interaction  in  the  nozzle  collar.  Thrust 
efficiencies  indicated  that  whistling  could  be 
achieved  without  much  greater  losses  than  the  basic 
axisynmetric  jet.  Entrainment  tests  were  inconclusive 
regarding  the  whistler  nozzle  performance.  Second, 
supersonic  flow  past  an  oscillating  cylindrical  shell 
is  analyzed  using  linearized  characteristics  methods. 
Pressure  distributions  and  generalized  aerodynamic 
forces  are  calculated  and  presented  for  various  radius 
to  length  ratios  and  reduced  frequencies.  Good 
agreement  is  ohtained  in  the  two-dimensional  limiting 
case  with  previous  work  by  Platzer,  and  an  early 
solution  of  the  steady  cylindrical  case  by  Zierep. 
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I-   IOEQDUCTI0N-THE  WHISTLER  NOZZLE 

The  phenomenon  of  entrainment  is  a  process  through  which  a 
moving  fluid  increases  its  mass  flow,  by  "grabbing  up"  the 
medium  it  is  flowing  in.  The  moving  fluid  is  called  the 
working  fluid,  and  might  be  the  primary  water  jet  in  an 
ejector  pump.  The  medium  might  be  any  other  fluid,  say  oil, 
in  the  case  of  the  ejector  pump. 


The   entrainment   process 
quite   a   long   time 
entrainment  process  exi 


has  found  practical  use  for 
Analytical  descriptions  of  the 
t:  however, none  is  general  enough  to 

E v e n    the 


describe   all   cases   of   practical   interest. 


steady   jet  in  still  air  is 

which 


relatively   simple   case   of   a 

rather  difficult  because  of  the   many   configurations 

uid  like  to  describe.  Thus  interest  and  research  in 
continues,  perhaps  with  more  intensity  than 
previously  due  to  the  recent  interest  in  new  V/STOL 
configurations  which  would  employ  this  process.  One  such  is 
Havy's   XFV  -  12A,   which   is  expected  to  derive  a  large 

or   entrainment, 


one   wo 
this   area 


the 


amount  of  thrust  through  the  augmentation, 

of  its  primary  thrust  using  wing  mounted  ejectors, 


The  success  of  concepts  such  as  the  augmenter  wing 
depends  en  many  factors,  an  important  one  being  the  rate  of 
entrainment  cf  the  primary  jet,  and  another  being  its 
efficiency.  This  rate  of  entrainment  is  the  rate  of 
increase  of  jet  mass  flow  with  axial  distance  away  from  the 
jet  nozzle.  The  rate  of  entrainment  and  nozzle  efficiency 
are  characteristics  which  vary  with  nozzle  configuration. 


As   mentioned   earlier   ,   there   is   a  large  variety  of 
steady  nozzle  configurations.   Lobe,  slot,   and   hypermixing 


nozzles  are  some  which  have  drawn  recent  interest  for  their 
entrainraent  capabilities.  One  should  not  neglect  unsteady 
nozzle  configurations,  however,  as  such  effects  as  swirling 
the  primary  flow,  and  oscillating  the  primary  flow  in  a  bi- 
stable manner  (fluidic  oscillator),  also  show  promise  in 
their  ability  to  entrain. 

One  of  the  newest  nozzle  configurations  to  surface  is  an 
unsteady  nozzle.  Perhaps  the  simplest  in  construction  of 
all  unsteady  nozzles,  it  is  illustrated  in  figure  1.  The 
whistler  nozzle,  as  it  is  called,  produces  a  loud  pure  tone 
while  in  operation.  Additionally,  it  creates  a  relatively 
large  rate  of  entrainment.  It  is  the  whistler  nozzle  which 
is  the  subject  of  this  thesis. 
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Figure    1  -   WHISTLER  NOZZLE  CROSS  SECTION 
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II.   RELATED  PAST  WORK 

The  earliest  work  done  on  the  v;histler  nozzle  was  by  Hill 
and  Greene  (reference  1) .  In  1974,  they  cited  the  discovery 
of  a  new  phenomenon,  characteristic  of  a  new  device  which 
they  named  the  whistler  nozzle.  Their  basic  configuration, 
as  illustrated  in  figure  1,  consisted  of  a  convergent 
section,  followed  by  a  section  of  constant  area  followed  by 
a  step  change  in  cross  sectional  area  to  another  section  of 
constant  area.  The  length  of  the  last  section  could  be 
varied  tc  produce  loud  pure  tones  of  varying  frequency, 
through  the  excitation  of  a  standing  acoustic  wave  in  the 
final  section,  much  as  in  an  organ  pipe.  Hill  and  Greene- 
observed  that  the  whistler  produced  a  greatly  increased 
mixing  rate,  and  decided  that  this  was  due  to  acoustic 
stimulation  of  the  jet.  Increased  mixing  rate  of  an 
axisymmetric  jet  due  to  acoustic  stimulation  is  a  process 
that  had  been  observed  earlier  by  Crow  and 
Champagne  (reference  5). 

In  addition  to  the  basic  axisymmetric  configuration, 
many  other  configurations  were  found  to  produce  the 
whistling  phenomenon.  A  few  are  illustrated  in  figure  2. 
Common  to  all  of  these  was  that  the  step  in  cross  section 
had  to  extend  completely  around  the  jet.  Also  noted  was  the 
apparent  presence  of  a  separation  reattachment  cycle, 
occurring  at  the  nozzle  lip. 

Kc  further  study  seems  to  have  been  performed  concerning 
the  whistler  nozzle,  save  a  report  in  1975  by  Hill  and 
Jenkins  (reference  2) ,  in  which  they  reestablished  the 
operating  characteristics  of  the  nozzle,  and  its  ability  to 
entrain  more  air. 
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Hill  and  Greene's  work  has  been  summarized,  since  their 
original  report,  in  comprehensive  reviews  of  the  ejector 
state  of  the  art  by  Schura  (reference  3) ,  and  by  Viets 
(reference  4).  These  reports  call  attention  to  the  whistler 
nozzle's  promise  in  rate  of  entrainment,  and  the  lack  of 
data  regarding  the  whistler  thrust  efficiency. 

Although  no  further  research  directly  involving  the 
whistler  is  available,  much  work  has  been  done  in  the  area 
of  free  jets,  and  entrainment,  which  might  prove  valuable  in 
a  study  of  the  whistler  nozzle. 

An  obviously  important  area  of  related  work  is  the 
sensitivity  cf  free  jets  to  acoustic  stimulation.  PrcbaDly 
the  best  known  and  most  definitive  work  in  this  area  was 
performed  by  Crow  and  Champagne  (reference  5) .  By  mounting 
a  speaker  inside  the  plenum  chamber  they  managed  to  excite  a 
free  axisymmetric  jet,  at  various  frequencies.  The  effect 
of  the  sound  was  to  produce  a  velocity  fluctuation  in  "cne 
nozzle  exit  plane.  It  was  observed  that  the  turbulent 
structure  of  the  jet  was  sensitive  to  the  acoustic 
excitation.  Through  the  use  of  high  —  speed  schlieren 
photography  it  was  observed  that  under  acoustic  stimulation 
the  turbulent  structure  of  the  jet  developed  a  large-scale 
vortex  structure,  as  contrasted  with  the  small  scale  vortex 
sheet  surrounding  a  normal  jet.  Attendant  to  the  formation 
of  the  large-scale  vortices,  an  increase  in  the  mixing  rate 
was  also  observed. 

Others,  such  as  Bevilaqua  and  Lykoudis  (reference  6) , 
have  observed  the  turbulent  structure  in  steady,  free  jets. 
The  presence  of  a  large  -  scale  vortex  structure  has  been 
observed  in  such  jets  at  low  Reynolds  numbers.  It  appears 
that  the  large  —  scale  structure  dissappears  as  Reynolds 
number  increases.   However,  according   to   the   findings   of 


14 


Crow  and  Champagne,  it  may  be  reexcited  through  accusiic 
stimulation.  It  may  be  conjectured  from  experimental 
evidence  that  the  entrainment  process  is  a  combination  of  a 
large—scale  "scooping  up"  of  the  ambient  fluid,  and  a  small- 
scale  nibblir.g  away  at  the  ambient  fluid.  These  entraiument 
processes  wculd  correspond  to  the  large  and  small  scale 
vortex  structures  around  a  free  jet. 

A  further  area  of  related  research,  perhaps  not  obvious, 
is  the  well  known  Coanda  effect.  The  basic  whistler 
configuration  sets  up  a  flow  situation  somewhat . aKin  tc  that 
of  a  bistable  fluidic  element,  or  fluidic  oscillator.  The 
difference  between  the  two  cases  is  that  the  fluidic 
oscillator  oscillates  between  two  positions,  while  the 
whistler  must  oscillate  between  many  circumferential 
positions.  That  such  an.  oscillation  takes  place  seems 
indicated  by  the  nozzle  configuration  and  the  separation 
reattachment  cycle  noted  by  Hill  and  Greene. 

Reference  10  provides  a  good  summary  of  the  mechanism  of 
the  Coanda  effect,  as  well  as  its  characteristics  regarding 
reattachment  and  sensitivity  to  sound.  Though  the  following 
comments  relate  specifically  to  two-  dimensional  flows,  it 
is  thought  that  the  general  behavior  and  mechanism  should 
carry  over  to  a  three-  dimensional  situation. 

Figure  3  illustrates  the  Coanda  mechanise  in 
two-dimensional  flow.  As  shown,  a  jet  is  issuing  frcn  an 
orifice  and  flowing  between  two  adjacent  walls.  Because  the 
jet  will  at  some  time  be  nearer  to  one  wall  than  the  ether, 
and  because  the  jet  is  entraining  air,  a'  pressure  difference 
across  the  jet  will  develop  that  pulls  the  jet  to  the  closer 
wall.  This  attachment  forms  a  region  of  circulating  flow. 
The  point  of  jet  attachment  is  determined  by  this  region,  in 
which  a  balance  is  struck  between  flow  entrained  by  the  jet 
and  flow  reinjected  by  the  jet. 
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Anything   that   disrupts   this  balance  causes  the  attachment 
point  of  the  jet  to  shift. 

Observation  has  shown  that  the  attachment  point  in 
Ccanda  effect  is  sensitive  to  jet  velocity  and  sound,  within 
a  range  of  Reynolds  numbers.  In  the  range  of  Reynolds 
number  from  200  to  3000,  under  the  effect  of  increasing  jet 
velocity,  or  a  sound  input  of  increasing  amplitude,  the 
attachment  point  has  been  observed  to  move  upstream.  This 
upstream  movement  appears  to  approach  a  limiting  position, 
beyond  which  the  point  will  not  move. 

Reference  10  relates  the  effects  of  jet  velocity  and 
sound  on  ih e  attachment  point  to  the  point  of  transition 
from  laminar  to  turbulent  flow  of  the  jet.  Both  increasing 
jet  velocity  and  a  sound  input  will  move  the  transi-ion 
pcint  closer  to  the  jet  exit.  This  movement  will  increase 
the  length  of  turbulent  flow  across  the  recirculation 
region.  A  turbulent  jet  entrains  more  air  than  a  laminar 
jet,  and  thus  the  balance  of  flow  in  the  recirculation 
region  is  disrupted,  with  a  resulting  upstream  move  of  the 
attachment  pcint. 
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III.   SCOPS  OF  RESEARCH 


The  intent  of  this  research  was  to  provide  a  clearer 
insight  into  the  workings  of  the  whistler  nozzle,  and  to 
develop  an  improved  understanding  of  the  ways  in  which  one 
may  effect  control  over  the  characteristics  of  a  free  jet. 

Specific  areas  of  research  were  the  nozzle  thrust 
efficiency,  the  nozzle  entrainment  rate,  and  visualization 
of  the  lip  interaction  process  described  by  Kill  and  Greene. 
Accordingly,  experiments  were  conducted  to  obtain  relevant 
data.  Thrusts  were  measured  on  a  thrust  bed,  mass 
entrainment  measured  in  a  device  patterned  after  that 
developed  by  Kicou  and  Spalding  (reference  7)  ,  and  high 
speed  motion  pictures  taken  of  tufts  about  the  nozzle  lip. 
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IV. 


NOZZLE  DESCRIPTION 


The  whistler  nozzle  used  in  all  tests  is  illustrated  in 
figure  4.  It  consisted  of  a  small  plenum,  convergent 
section,  short  constant  area  section  and  a  moveable  collar 
over  the  constant  area  section^  The  inner  whistler  diameter 
was  1  inch,  while  the  collar  diameter  was  1.5  inches.  The 
nozzle  configuration  was  essentially  that  used  by  Hill  and 
Greene  in  their  investigations.  The  characteristic 
dimension  used  for  Reynolds  number  and  nondimensionalization 
of  axial  distance  from  the  nozzle  was  the  inner  nozzle 
diameter,  1  inch. 


19 


Whistler 
Nozzle 


Plenum 


Presure  Tap 


Air  Inlet 


Figure    U  -   WHISTLER  NOZZLE  AND  PLENUM 
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v-   DESCRIPTION  OF  EXPERIMENTS 


A.   THRUST  EFFICIENCIES 


Thrust  efficiencies  were  obtained  by  measuring  nozzle 
thrust.  A  comparison  was  made  between  measured  and 
calculated  thrust. 

Test  equipment  consisted  of  a  thrust  bed  and  the  nczzle. 
The  thrust  bed  was  configured  to  measure  horizontal  thrust, 
and  slid  on  two  tracks,  mounted  to  a  solid  foundation.  The 
nozzle  plenum  was  clamped  to  the  thrust  bed  which  was  in 
turn  secured  to  the  foundation  through  a  lead  ring,  mounted 
with  a  load  cell,  the  device  that  actually  measured  thrust. 
Air  was  supplied  to  the  nozzle  from  a  large  tank,  connected 
to  the  plenuir  through  a  regulator  and  a  flexible  tube  which 
led  vertically  away  from  the  plenum,  so  as  not  to  bias  the 
thrust  measurements  in  any  way.  Calibration  was  performed 
with  the  use  of  a  tray,  connected  to  the  bed  through  a  cable 
and  pulley  arrangement,  upon  which  one  placed  weights  of 
known  magnitude.  Peripheral  gear  required  included  a 
voltage  source,  digital  voltage  meter,  and  Wheatstone 
bridge.  Figure  5  illustrates  the  set  up.  Plenum  pressures 
were  measured  using  a  mercury  manometer. 


21 


Figure    5    -       TFRflST    PED 
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Figure  6  -   THRUST  BED  CRLTEFSTIOK  CURVE 
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Calibration  of  the  thrust  bed  indicated  a  certain  amount 
of  hysteresis  in  the  system,  as  may  be  seen  in  figure  6. 
Not  serious  at  thrusts  of  about  6  pounds,  it  did  introduce 
problems  at  low  thrust  levels,  and  precluded  reliable  thrust 
measurement  at  low  Reynolds  numbers. 

Thrust  measurements  were  made  for  various  collar 
positions,  over  a  range  of  nozzle  exit  velocities.  Recorded 
data  included  gage  pressure  in  the  nozzle  plenum,  thrust 
readings  from  the  voltmeter,  collar  extension,  and  ambient 
pressure  and  temperature. 

Measured  thrusts  were  plotted  against  pressure  ratio 
across  the  nozzle.  Figure  7  presents  an  example  curve. 
Ideal  thrust  was  calculated  according  to  the  procedures  in 
appendix  A.  Thrust  efficiencies  were  then  obtained  and 
plotted  against  nozzle  pressure  ratio,  exit  velocity,  and 
Reynolds  number  (based  on  the  nozzle  diameter) . 


B. 


MASS  EN1EAINM2NT 


Mass  entrainment  was  measured  directly,  using  a  device 
patterned  after  one  developed  by  Ricou  and  Spalding 
(reference  7).  It  was  believed  that  this  device  should 
yield  more  accurate  results  than  an  approach  using 
integrated  velocity  profiles,  especially  in  the  case  of  an 
unsteady  nozzle.  Measurements  were  made  at  various  axial 
positions,  for  several  collar  positions. 

The  device  used  will  hereafter  be   referred  to   as   the 

entrainment   chamber.    It  was  designed  with  the  aid  of  data 

from  a  research  report  by  Peschke  (reference  8) .  Appendix  B 
provides  a  detailed  account  of  the  design. 
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figure  7  -   EXPERIRENTAL  THRUST  CURVE 
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The  principle  cf  operation  is  relatively  simple.  The 
entrained  air  of  a  free  jet  flows  in  at  a  right  angle  to  the 
jet  axis.  The  entrainment  chamber  simply  blocks  this  inflow 
witn  a  circumferential  membrane,  and  seals  the  ends  with  a 
plate  at  the  nozzle  and  an  orifice  at  the  axial  station 
where  one  wishes  to  measure  the  total  mass  flow.  Ideally, 
all  of  the  entrained  air  is  supplied  to  the  exterior  of  the 
membrane  from  a  regulated  source.  As  a  free  jet  does  not 
support  any  axial  pressure  yradient,  the  proper  amount  of 
air  to  supply  through  the  membrane  is  indicated  by  the 
disappearance  of  any  pressure  gradient  across  the  orifice. 
One  may  regulate  the  membrane  flow  to  obtain  this  point. 
These  primary  aspects  of  the  device  are  illustrated  in 
figure  8. 

Critical  design  aspects  of  the  entrainment  chamber, 
especially  the  orifice  sizes,  are  treated  in  appendix  B.  It 
is  sufficient  here  to  simply  state  that  one  entrainment 
chamber  was  designed  to  perform  measurement  of  the  total  jet 
mass  flow  at  various  axial  locations. 

One  primary  difference  between  this  chamber's  geometry 
and  that  of  ethers  (references  7  and  8)  should  be  noted.  As 
shown  in  figure  8,  exit  orifices  were  positioned  in  the 
porous  cylinder  at  various  distances  from  the  nozzle.  That 
portion  cf  the  porous  cylinder  extending  above  the  orifice 
plate  was  sealed  by  a  fiberglass  cylinder  attached  to  the 
orifice.  Thus, in  this  situation  that  portion  of  the  jet 
beyond  the  orifice  and  inside  the  porous  cylinder  exit  was 
not  in  a  truly  free-jet  flow  situation.  This  was  as  opposed 
to  the  reference  7  and  8  configurations  where  the  exit 
orifice  corresponded  to  the  exit  of  the  porous  cylinder. 
Thus,  in  those  experiments  the  jet  past  the  orifice  was  in  a 
truly  free-jet  flow  situation.  The  configuration  used  in 
this  report  was  adopted  in  order  to  facilitate  the 
measurement  of  mass  flows  at  many  axial  locations.    It   was 
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hoped  that  this  chamber's  geometry  would  not  affect  the 
chamber  readings  greatly,  although  some  effect  was  to  be 
expected . 

Peripheral  equipment  to  the  basic  entrainrnent  chamber 
included  flew  regulating,  metering,  aad  measuring  devices. 
Figure  9  presents  a  photograph  of  the  chamber  and  its 
peripheral  gear.  Separate  air  sources,  both  fed  by  a  large 
tank,  were  supplied  for  the  nozzle  and  membrane.  Nozzle 
flow  was  metered  through  an  orifice  and  controlled  with  a 
regulator.  Eecause  of  the  much  greater  mass  flow  required, 
membrane  air  was  supplied  through  a  venturi,  followed  again 
by  a  flow  regulator.  Pressures  across  the  oririoe  were 
measured  with  an  inclined  water  manometer.  Pressures  across 
the  venturi  were  measured  with  a  mercury  manometer. 
Pressure  across  the  entrainrnent  chamber  orifices  was 
measured  v,ith  an  alchohoi  micromanometer.  Nozzle  plenum 
pressure  was  measured  with  a  water  manometer. 
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Figure    9    -      THE    ENCHAINMENT    CHAM3EB    AND    PERTPHERAL    GEA:? 
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Membrane   air   temperature   was   measured   just   after    the 
venturi,  using  a  mercury  thermometer. 

Entrainment  measurements  were  made  at  various  axial 
locations  for  various  collar  positions.  Recorded  data 
included  pressure  difference  across  the  venturi  and  orifice, 
pressure  difference  across  the  entrainment  chamoer  orifice, 
nozzle  plenum  pressure,  ambient  pressure,  air  temperature  in 
the  membrane  supply,  and  collar  position.  Standard 
procedure  for  eacn  run  was  to  vary  the  membrane  airflow  so 
that  one  oftained  both  positive  and  negative  pressure 
gradients  at  the  chamber  orifice.  This  provided  a  more 
precise  location  of  tne  proper  entrained  flow  rate.  It  is 
of  interest  to  note  that,  due  to  hardware  limitations  and 
the  large  ancunt  of  air  entrained  by  the  jet  at  the  furthest 
axial  locations  measured,  nozzle  exit  velocities  were 
limited  to  around  100  feet  per  second  at  the  further  axial 
locations.  This  also  affected  the  number  of  useful 
locations  for  whistler  measurements,  as  a  strong  whistler 
oscillation  required  a  relatively  large  exit  velocity,  say 
300  feet  per  second. 

Data  reduction  followed  standard  procedures  for  the 
venturi  and  orifice.  These  are  presented  in  appendix  C. 
Entrainment  was  plotted  as  M/MO (total  mass  flow  divided  ny 
the  primary  mass  flow) ,  against  X/D,  axial  location  divided 
by  nozzle  diameter.  Working  curves  of  chamber  orifice 
pressure  gradient  are  presented  as  pressure  difference 
against  H/::0,  in  the  example  plotted  in  figure  10. 
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c, 


FLOW  VISUALIZATION 


Flow  visualization  of  the  nozzle  jet  lip  interaction  was 
performed  in  two  ways:  high-speed  motion  photography  of  a 
ring  of  tufts  at  the  nozzle  lip,  for  the  purpose  of 
discerning  a  pattern  to  the  interactions;  and  a  mixture  of 
lamp  black  end  kerosene  inside  the  nozzle  collar  in  order  to 
observe  the  behavior  of  the  separation  reattachment 
location . 

Photographic  equipment  consisted  of  a  Honeywell  Fentax 
single  reflex  camera,  a  Eolex  Supreme  handheld  16mm  movie 
camera,  and  a  ftedlake  Hycam  16mm  high-speed  movie  camera. 
The  tuft  ring  was  simply  constructed  of  cotton  threads  held 
by  a  masking  tape  ring,  and  a  photo  of  the  ring  is  presented 
in  figure  11.  Observation  in  the  lampblack  experiment  was 
performed  with  the  naked  eye  and  no  peripheral  equipment  was 
necessary. 

Photographic  tests  were  to  be  done  in  three  phases. 
High-speed  still  shots  were  first  taken  to  determine  the 
usefulness  of  the  tufts  for  visualization,  and  to  obtain  the 
approximate  required  film  speed  for  motion  photography. 
Next,  mediutD  speed  motion  pictures,  64  frames  per  second, 
were  taken  to  confirm  the  necessity  of  the  high-speed 
camera.  Finally,  high-speed  motion  pictures  were  taken  to 
visually  control  the  tuft  motion  such  that  useful 
observations  could  be  made.  No  particular  collar  settings 
or  exit  velocities  were  used  in  photographic  work.  The  only 
requirement  was  to  produce  the  oscillation  phenomenon,  which 
was  never  difficult.  As  it  turned  out,  a  lack  of  time 
prevented  the  completion  of  the  high  —  speed  photography 
phase. 
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Lacipblack  tests  were  performed  in  two  ways.  The  inside 
of  the  cciiar  extension  was  liberally  coated  with  the 
lampblack  mixture,  ar.d  either  exit  velocity  varied  at 
constant  cellar  setting,  or  collar  setting  varied  at 
constant  exit  velocity.  Observations  were  all  visual,  none 
quantitative,  save  noting  the  range  of  exit  velocities 
covered,  which  was  from  about  100  to  300  feet  per  second. 
The  lampblack  mixture  was  observed  to  pool  as  illustrated  in 
figure  12,  and  the  downstream  —  most  edge  of  the  pool  was 
interpreted  as  indicating  the  location  on  the  cellar  of 
reattachment  and  separation. 
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Fiaure  12  -   LAMPBLACK  POOLING  ILLUSTRATION 
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VI.   RESULTS  AND  DISCUSSION 


Nozzle  thrust  efficiency  results  are  presented  in  figure 
13,  plotted  for  various  collar  positions.  The  whistler 
oscillation  was  produced  at  each  collar  setting  for  the 
whole  range  of  Reynolds  numbers  shown,  although  sudden 
frequency  jumps  did  occur  in  this  range.  These  frequency 
jumps  did  net  appear  to  change  the  nozzle  thrust  or  thrust 
efficiency  in  any  discontinuous  manner.  Instead  the  trend 
was  as  indicated,  progressively  decreasing  nozzle  efficiency 
with  increasing  collar  extension.  The  zero  collar  extension 
efficiencies  were  found  to  oe  about  6  percent  lower  than 
those  found  in  a  previous  report  (reference  8)  .  This  is 
believed  due  to  the  extremely  simple,  unflared  configuration 
of  the  nozzle,  which  may  be  seen  in  figure  U.  Reference  8 
nozzles  were  divergent  nozzles  and  as  such  should  be  more 
efficient.  Important  to  note  is  that  the  whistler 
oscillations  may  be  obtained  without  much  reduction  in  the 
basic  nozzle  efficiency. 

Entrainment  results  are  presented  in  figure  14,  plotted 
together  with  the  results  of  Hill  and  Greene  (reference  1) 
for  comparison.  As  may  be  seen,  in  the  case  of  the  basic 
jet  agreement  is  good  qualitatively  but  poor 
quantitatively.  In  the  case  of  the  whistler  nozzle 
agreement  is  not  good  in  any  way.  It  is  believed  that  the 
poor  quantitative  comparison  was  due  to  interference  of  the 
deep  orifice  cylinders  with  the  jet.  These  create  a 
pressure  gradient  from  the  orifice  to  the  chamber  exit  which 
quite  possibly  affects  the  readings.  As  noted  in  the 
experiment  section,  and  in  Appendix  C,  these  orifice 
cylinders  constitute  a   major   difference   in   configuration 
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from  the  chambers  of  reference  7  and  8,  Additionally,  in 
the  case  cf  the  whistler  nozzle,  it  is  believed  the  method 
of  enclosing  the  jet  up  to  the  axial  position  of  measurement 
may  interfere  with  the  sound  -  jet  interaction  (especially 
the  cloth  porous  cylinder)  and  the  formation  of  a  larcj^- 
scale  turbulent  structure.  Figure  15  illustrates  that 
readings  taken  were  relatively  constant  with  Reynolds 
number.  As  noted  by  Ricou  and  Spalding,  the  entrainment 
ratio  can  be  a  strcng  function  of  Reynolds  number  in  certain 
ranges  of  Reynolds  number.  It  is  interesting  to  note  that 
this  behavior  is  supported  by  the  present  notion  of  jet 
turbulence  structure  and  its  development  as  exit  velocity 
increases,  which  postulates  that  the  major  contributor  to 
entrainment,  the  large  scale  structure,  dissappears  as 
Reynolds  numrer  increases. 
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Ficrure       13    -       THRUST    EFFICIENCY    RESULTS 
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Figure   14  -   ENTRAPMENT  RATIO  RESULTS  AND  COMPARISON  DATA 
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Lack  of  time  precluded  testing  at  more  jet  axial 
locations,  as  each  location  required  a  different  orifice 
cylinder  (see  appendix  A).  Of  the  original  four  orifice 
cylinders  constructed  using  reference  8  data,  only  two  were 
useable  as  the  spreading  angles  encountered  averaged  about 
28  degrees,  cr  10  degrees  less  than  those  used  in  reference 
8.  In  order  to  be  effective,  the  orifices  had  to  be 
positioned  at  axial  distances  further  than  they  were 
designed  fcr.  One  of  those  was  modified  to  obtain  readii 
at  its  design  position,  also .  The  larger  two  were  never- 
effective.  Another  problem  encountered  was  the  mass  flow 
limitation  cf  the  porous  cylinder  air  supply.  This 
restricted  the  useable  nozzle  flow  at  axial  locations  beyond 
X/B  of  about  6  to  such  low  values  that  a  whistler 
oscillation  strong  enough  to  affect  the  entrainment  could 
not  be  produced. 

It  was  hoped  that  high  speed  notion  pictures  i.ould  b 
taken  of  the  tuft  experiment;  however,  this  phase  of  the 
photographic  experiments  was  not  reached.  The  effectiveness 
of  the  tufts  in  indicating  the  circumferential  collar 
positions  of  attached  and  detached  flow  was  proven  in  the 
high-speed  still  picture  phase.  Motion  pictures  with  a  b^i 
frame  per  second  camera  showed  that  a  much  higher  speed 
camera  was  necessary  to  slow  the  motion  such  that  motion 
pictures  would  be  useful  in  slowing  the  tuft  motion  to  a 
reasonable  speed.  The  necessary  film  speed  was  estimated  as 
about  800  f rarces  per  second;  however,  the  pictures  were  never 
taken.  Seme  representative  still  pictures  are  presented  in 
figure  16. 

Lampblack  tests  appeared  to  be  productive  as   indicators 
of   the   general   behavior   of   the  collar  attachment  point. 
Indications  are  that  a  sort  of  Coanda   effect   is   involvedj 
modified   by   the  fact  that  while  the  jet  is  attached  at  any 
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point  on  the  collar  there  are  other  points  where  detach 
flow  allows  a  backflow  into  the  collar  jet  gap, which  tends 
to  relieve  the  Coanda  effect  causing  a  circumferential 
shift  in  the  point  (or  points)  of  attachment .  At  a  fixed 
collar  extension,  with  increasing  jet  velocity,  the  attachment 
pcint  was  observed  to  move  upstream  to  a  limiting  position, 
just  as  the  Coanda  jet  described  in  reference  10.  However, 
at  fixed  exit  velocity  and  increasing  collar  extension  a 
different  behavior  was  observed.  The  attachment  point 
seemed  to  move  with  the  collar  until  a  point  was  reached  at 
which  the  frequency  of  oscillation  jumped.  At  that  instant 
the  point  of  attachment  also  appeared  to  jump  suddenly  back 
to  its  original  position  relative  to  the  jet  exit.  This 
behavior  indicated  some  interdependence  of  the  v.'histler 
attachment  phenomenon  and  the  oscillation  phenomenon. 

In  addition  to  the  above  experiments,  two  other 
observations  were  made  which  verified  those  of  Hill  and 
Greene. 

First,  it  was  noted  that  extension  of  the  whistler 
collar  seemed  to  increase  mass  flow  through  the  nozzle. 
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Figure   16  -   TUFT  PICTURES,  RESULTS 
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This  was  confirmed  with  static  pressure  measurements  in  the 
nozzle  'throat1.  These  indicated  that  tne  collar  was  acting 
as  an  overexpanded  diffuser,  resulting  in  an  effective  exit 
area  slightly  greater  than  the  'throat'  cross  section  but 
less  than  the  collar  cross  section,  and  further  resulting  in 
a  greater  mass  flow  through  the  nozzle. 

Second,  it  was  noted  that  a  correspondence  between 
frequency  of  oscillation,  jet  exit  velocity,  and  collar 
extension  seemed  to  exist.  It  was  found  that  by  gradually 
extending  the  collar  as  plenum  pressure  was  increased,  one 
could  'follow*  a  frequency  of  oscillation.  Accordingly, 
collar  extensions  were  measured  and  jet  velocity  calculated 
at  fixed  frequencies  of  oscillation,  and  it  was  found  that 
these  fixed  frequencies  corresponded  to  constant  ratios  of 
collar  extension  to  jet  velocity,  thus  confirming  the  organ 
pipe  character  of  the  whistler  oscillation  noted  by  Kill  and 
Gre  ene. 
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VII.   CONCLUSIONS  AND  RECOMMENDATIONS 


The  thrust  efficiency  and  mass  entrainment 
characteristics  of  the  whistler  nozzle  have  been  measured, 
in  the  case  cf  the  entrainment,  using  a  technique  heretofore 
not  used  on  an  unsteady  jet.  Further,  observations  have 
been  made  which  should  contribute  to  a  further  understanding 
of  the  whistler  flew  mechanism. 

Several  directions  of  further  work  are  indicated. 
First,  it  is  recommended  that  a  thrust  measurement  device 
using  a  pendulum  arrangement,  similar  to  that  used  in 
reference  15f  be  constructed  for  a  more  accurate 
determination  of  the  whistler  efficiency.  A  device  of  this 
type  might  avoid  the  hysteresis  problems  encountered  with 
low  thrusts  on  a  conventional  thrust  bed.  It  is  also 
recommended  tnat  a  whistler  nozzle  of  smaller  throat 
diameter  be  constructed  to  permit  mass  entrainment 
measurements  ,  at  further  jet  axial  locations,  in  the 
entrainment  chamber,  at  the  high  exit  velocities  required  to 
produce  a  strong  oscillation.  This  will  also  permit  further 
study  of  the  possibility  of  whistler  measurements  in  the 
chamber  ever  a  greater  range  of  Reynolds  number.  Next  it  is 
recommended  that  further  flow  visualization  experiments  be 
conducted.  Specifically,  high  speed  motion  pictures  of  the 
tufts  used  here,  and  high  speed  schlieren  photography  cf  the 
jet,  should  prove  useful.  Finally,  it  is  recommended  that 
other  unsteady  nozzles  be  tested  in  the  chamber,  as  the  poor 
results  obtained  here  with  the  whistler  are  only  preliminary 
and  may  not  te  indicative  of  the  usefulness  of  the  chamber 
with  non-acoustic  effects  such  as  swirl. 
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VIII.   INTRODUCTION-  THE  RING  WING 


The  present  work  concerns  a  hollow  cylinder,  axially 
aligned,  in  supersonic  flow  and  undergoing  small  amplitude 
oscillations  in  angle  of  attack.  The  intent  is  to  calculate 
the  resulting  pressure  distributions  and  aerodynamic  moments 
on  what,  sight  be  mere  properly  termed  a  ring  wing. 

Previous  investigations,  i.e.  references  12  and  13,  have 
applied  the  method  of  characteristics  towards  the  separate 
calculation  or  the  inner  and  outer  flowfields,  as  well  as 
pressure  distributions  and  generalized  aerodynamic  forces, 
for  such  a  cylinder  in  supersonic  fiov:  whose  walls  are 
undergoing  small  amplitude  sinusoidal  oscillations,  or  panel 
flutter.  These  previous  investigations  resulted  in  computer 
programs  to  perforin  the  necessary  calculations. 

The  primary  thrust  of  the  present  work  was  the 
modification  of  the  previously  created  programs  and  then  the 
combination  of  their  results  to  calculate  the  unsteady 
pressure  distribution,  lift  and  moment  on  the  cylinder 
undergoing  angle  of  attack  oscillations.  Finally,  using  the 
modified  programs,  a  parametric  study  of  these  items  was 
conducted. 
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IX- 


PROBLEM  FOR^uLnriON 


This  development  is  closely  based  on  several  previous 
papers  ,  i.e.  references  12  and  13. 

Consider  a  circular  cylindrical  shell,  in  supersonic, 
inviscid,  adiabatic  flow,  whose  axis  is  aligned  with  the 
freestreara,  and  whici:  is  undergoing  small  amplitude 
oscillations  in  angle  or  attack,  about  the  Z  axis,  (rigure 
17) 

The  governing  equation  for  this  flow  situation  is  the 
linearized   unsteady   potential   equation,   in   cylindrical 

coordinates. 


(1-M2)  $    +  $»    +  1  $   +1   $, 
v     '   xx     rr       r     2 

r      r 


2M  $. 


1   »    =  0 

„2    tt 


The  relevant  boundary  conditions  are  the  flow  tangency 
condition  and  Sommerf eld1 s  radiation  condition,  i.e., that 
disturbances  will  propagate  away  from  their  source. 


m 


X 


+J 

rH 


A 


Figure   17 


RING  WING  ILLUSTRATION 
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One  further  condition  is  necessary,  along  the  X-axis,  which 
amounts  to  setting  the  pressure  disturbance  at  the  axis 
equal  to  zero. 

The  linearized  form  of  the  flow  tangency  condition  is: 


'rl     =  lh  +  U   8h  ,0  N<  x  $  l 
r=R    3t     »  9x 


which  must  be  satisfied  on   the   cylinder   inner   and   outer 
surface. 


As  stated,  the  boundary  condition  at  the  axis  i 
pressure  disturbance: 


s  a   zero 


'P 


=   0 


r=0 


The  linearized  form  of  the  pressure   coefficient   is   as 
fellows: 


u^  t   Uro       qoo 


After  nondi&ensionalization,  using  cylinder  length  and 
length  divided  by  freestream  velocity,  the  above  eguations 
become: 
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(1-M2)     A  +    (J)  +14>       +    1 

XX  rr  —         r  —2 


OG 


"    2M2    *xt    "    M"    *tt    =    ° 


C_   =   -2(    di      +    <b      ) 

P  t  x 


subject    to: 


r^R         St         3x 


,0     i'     x     v<     1 


One  row  assumes  a  cylinder  oscillation  of  the  form 


h(x ,6 ,t)  =  Z (x)  cos (6)  e 


ikt 


where    the   non-dimensional    frequency,    kr    is: 


k  =  co_l 

Uco 


The  deflection   amplitude,   Z  (x) ,   for   the   case   under 
consideration,  is  of  the  form: 


Z  (x)  =  a  (  x  -  x   ) 

o 


50 


Following      the      assumed      form      of      the      deflection,      the 
perturbation    potential    must    be   of   the    form: 


(J)(x,r,0,t)    =    <t>  (x  ,r)    cos  (9)    e 


ikt 


and,  after  substitution,  the  basic  equations  become: 


(1-M2)    <fc        +    4>        +    1    <f)      +    (k2M2-l    )    <j>    -    i2kM   <})      =    0 
xx  rr         —      r ^r7- 


X 


C      =    -2(    (j>      +    ik6    ) 
P  x 


subiect    to: 


»  =    Z      +   ikZ         ,o   ^    x   ^    l 

r '  x 

r=R 


and : 


C    |=         0 
P   r=o 
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X.   «3TK0D  OF  CHARACTERISTICS 


For  supersonic  flow,  the  linearized  perturbation 
potential  equation  is  hyperbolic,  and  has  characteristics 
which  satisfy  the  following  differential  equation: 


(M2-l)  d2r  -  d2x  =  0 


Since  the  characteristics  in  supersonic  flow  correspond 
to  the  Mach  lines,  then  for  ds,  the  differential  arc  length 
along  a  characteristic: 


dr  =  +  _1  ,      dx  =  /m2-1 

ds      M  ds       M 


An    arbitrary    function,   F(x,r),   has   the   following 
derivatives  along  the  characteristics: 


dF  =  F   dx  +  F   dr  =  /Mz-1  Fv  +  1  F 
. ,     x  r  x  —  —  j^ 

ds       ds       ds       M       M 


Defining   ds1  and  ds2  to  be  the  differential  arc  lengths 

along  the  left,  and  right-running   Mach   lines   respectively, 

then: 

F1   =   dF       ,   F2  =  dF 
.ds-,  ds„ 
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The  cross  derivatives  become: 


12     "21 


F    =  1  (  (1-M2)  F    +  F   ) 
"^p         xx    rr' 


Solving    now   for    Fx   and   Fr: 


F.    =         M  (F      +    F.) 

x         2/M^l         1  2 


F      =   M     (F.     -    F~) 

r      T      1         2 


Now  the  basic  equations  may  be  written  in  terms   of 
and  its  derivatives  along  the  characteristics  as  follows: 


12 


=   1   (<f>.  -  (f)0)  +  (k: 
21    2rM    1     2 


1   )  (j>  -  ikM  (<f>   +  0  ) 
M2-l   -l     2 


r2M2 


<J)r  =  M  (<j)1  -  (|>2)  =  Zx  +  ikZ 


C   =  -2(4>   +  ik<|>) 
P        x 
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These  equations  have  been  put  into  finite  difference 
form  and  programmed  in  Fortran  IV  for  the  separate  cases  of 
exterior  flow  and  interior  flow  of  a  cylinder  whose  walls 
were  undergoing  small  amplitude  sinusoidal  oscillations.  In 
the  interior  flow  program  the  X  axis  boundary  condition  is 
satisfied  on  a  small  inner  cylinder  in  place  of  the  axis,  to 
circumvent  problems  encountered  should  r  go  to  zero.  For 
detailed  accounts  of  the  equations  and  programs  refer  to 
references  12  and  13. 

The  primary  work  here  was  the  transformation  of  the 
af orementioned  programs  to  the  case  under  consideration  , 
i.e.,  a  cylinder  undergoing  small  amplitude  oscillations  in 
angle  of  attack,  and  then  a  parametric  variation  study  of 
the  pressure  distributions  and  generalized  aerodynamic 
force.  The  programs  were  modified  such  that  all  cases  run 
were  with  an  angle  of  attack  amplitude  of  0.1  radians. 
Since  the  basic  eguation  is  in  a  linearized  form,  the 
conversion  of  any  results  to  other  angle  of  attack 
amplitudes  is  a  simple  matter  of  multiplication.  For  those 
interested,  the  programs  resulting  from  changes  required  to 
the  reference  12  and  13  programs  are  listed  in  appendix  D. 
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XI.   GENERALIZED  AERODYNAMIC  FORCE 


As  stated  above,  the  results  studied  were  the  pressure 
coefficient  amplitude  distribution,  and  the  generalized 
aerodynamic  force,  Qmr.  In  its  most  general  form,  Omr  is  as 
follows: 


Qmr  =  1_  /C    ¥r  (x)  dx 


Pm 


Z (x)  =  Z  A.  ¥ . (x) 
1       ~J       3 


where  Crm  is  the  pressure  coefficient  amplitude  resulting 
from  the  m'th  axial  deflection  mode,  PSIr  is  the  r'th  axial 
deflection  mode,  and  A  is  the  amplitude  of  the  j • th  axial 
deflection  mode.  For  the  special  case  under  consideration, 
Q  takes  the  following  f. or:  : 


Q  =  1  /C   ^ (x)  dx 
2        P 


Z  (x)  =  H»(x)  =  a(x  -  x  ) 
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Thus  ,  in  the  case  under  study  Q  has  the  special 
interpretation  of  being  the  moment  coefficient  amplitude 
about  the  point  Xc  multiplied  £sy  the  angle  of  attack 
amplitude,  alpha,  which  is  equal  to  0.1  radians  in  all  cases 
presented  . 
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XII.   RESULTS  AND  DISCUSSION 


Several  cases  cf  comparison  were  considered  towards  the 
evaluation  cf  program  results.  Figure  18  presents  two  cas 
of  comparison  for  pressure  distributions.  The  outer  radius- 
to-  length  ratio  was  taken  to  such  a  value  (10)  as  to 
approximate  a  flat  plate  locally.  The  inner 
radius~to-lengr.h  ratio  was  taken  to  such  a  value  (9.75)    as 


to 


c  a  u  s 


rpf 


reflection 


f  i 


on 


the  inner  surface  of  the  outer 


cylinder. 


For  sufficiently  large  radius  -  to  -  length  ratios  the 
exterior  pressure  distribution  should  approximate  the 
theoretical  value  for  a  fiat  plate  in  supersonic  flow,  i.e. 
2a/&.  Figure  18A  presents  this  comparison,  and  as  can  be 
seen   agreement  is  good. 

For  sufficiently  large  radius  -  to  -  length  ratios  the 
interior  pressure  distribution  should  approximate  the 
theoretical  value  for  a  flat  plate.  Additionally,  if  the 
inner  cylinder  is  of  such  a  radius-to-length  ratio  as  to 
cause  a  Mach  wave  to  reflect  upon  the  outer  cylinder  inner- 
surface,  a  flow  situation  exists  which  should  approximate  a 
flat  plate  in  a  free  jet,  where  the  mach  waves  are  reflected 
from  the  jet  boundary.  This  is  due  to  the  boundary 
condition  imposed  on  the  inner  cylinder,  Cp  =  0.  Reference 
16  has  treated  this  case  in  analytic  form  for  small 
frequencies.  Figure  18,  A  and  B,  presents  a  comparison  of. 
program  results  and  the  results  of  reference  16.  As  can 
seen  agreement  is  good. 
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Figure   18  -   AXIAL  PRESSURE  DISTRIBUTION,  2D  LIMITING  CASE 
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Reference  14  has  treated  the  cane  of  a  nonoscillati ng 
ring  wing  in  supersonic  flow.  Figure  19  presents  a 
comparison  of  procrram  results  and  those  of  reference  1'4,  for 
the  case  shown.   As  can  be  seen   agreement  is  fair. 

In  reference  12,  which  developed  the  basic  inner  flow 
field  program  used  in  this  report,  the  occurrence  of  a 
numerical  instability  in  the  pressure  distributions,  under 
certain  conditions,  was  noted.  it  was  found  that  this 
disturbance,  which  occurred  only  after  a  reflection,  war;  due 
to  an  undefined  constant  in  the  subroutine  which  corapn 
the  boundary  condition  at  the  axis.  After  defining  the 
constant  correctly,  the  pressure  distributions  were  se^n  to 
settle  down  cuitc  nicelv.  Figure  20  presents  a  before  Bnd 
after  comparison  of  pressure  distributions  for  a 
representative  case. 
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Figure   20  -   NUMERICAL  INSTABILITY  ILLUSTRATION 
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Figures  21  through  23  present  various  parametric  studies 
of  the  moirent  coefficient,  Q,  for  cases  Loth  with  and 
without  a  reflection  on  the  inner  cylinder  surface  In 
these  figures  the  coefficient  plotted  is  actually  the 
difference  letween  that  calculated  for  the  inner  cylinder 
surface  and  that  calculated  for  the  outer  cylinder  surface, 
divided  by  the  angle  of  attack  amplitude.  Thus  the 
coefficients  plotted  are  total  moment  and  pressure 
coefficients,  as  opposed  to  those  plotted  earlier,  which  are 
the  coefficients  for  the  inner  and  outer  cases  separately. 
As  can  be  seen  from  figure  22, for  a  center  of  rotation  at 
the  leading  edge  the  cylinder  has  definite  geometric 
regions  of  stability  with  respect  to  angle  of  attack. 
Figure  21  illustrates  the  effect  of  increased  frequency  on 
the  moment  coefficient.  Figure  23  illustrates  a  typical 
pressure  distribution  including  one  reflection,  from.  which 
one  can  infer  the  reason  for  the  behavior  of  Q  in  figure  22, 
and  also  the  behavior  of  Q  with  a  varying  center  of 
rotation. 
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XIII.   CONTIJJSIONS  AND  RECOMMENDATIONS 


The   computer  programs  of  references  12  aiul  13  have  been 
successfully   adapted   to   the   case   of   a   ring   wina    in 
supersonic   flow  undergoing  oscillations  in  angle  of  attack. 
Available  cases  for  comparison  have   shown   good   agreement, 
indicating   the   validity   of  the  method.   The  programs  !r 
possible  the  relatively  easy   calculation   of  the   unsteady 
lifts   and   moments   of   a  rinq  wing  in  supersonic:  flow,  for 
such   applications   as   tubular   artillery    shells    (STU0, 
spinning    tubular   projectile),  although   spin   nay   induce 
aerodynamic   effects   not   accounted   for   in   the    method. 
Further   evaluation   and   parametric  studies  are  recommen" 
for  a  more  complete  verification  o^  the  computer  programs. 
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APPTM-TPIX  A 


TSENTROPTC  THRUST  CALCULATIONS 


Isentropic   thrust   calculations   made   use   of  standard 
isentrooic  relations.   Calculations   wore   based   on   Die- 
pressure   and   nozzle  exit  area.   The  formula  for  isentropic 
thrust  is  arrived  at   3?   follows.    Nozzle   thrust   may   ^e 
written  in  the  following  form: 

T   =  w  V    =(pAM2/RT) 
s      ex  ex 

where  W  is  nozzle  mass  flow,  A  is  exit   area,   and   one   has 
applied   the   speed  of  sound  relation.   Exit  Hach  number  i 
be  written  in  terms  of  pressure  ratio  as  follows: 


y-i 


Vi    ,  =      2   ((P./P   )  '  V*  -l) 


ex 


jr-i 


ex 


Applying  this  and  the  state  equation  to  the  formula  for 
Ts  one  obtains: 

-1 

To  =  P^A     2  ((pr/p   >     -1) 
s     ex  ex  yzj        0   ex 

Now,  once  one  has  measured  plenum  and  ambient  pressure, 
and  thrust,  Texp,  one  can  calculate  the  isentropic  thrust 
efficiency: 


n  =  t   /  t 

exp'   s 
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APPENDIX  B 


ENTRAINMENT  RATIO  CALCnL.ITIONS 


Mass   flows   and   entrainment   ratios   were    calcula- 
according   to  the   equations   and   procedures   sat   out   in 
reference  11f  the  RSME  Power  Test  Codes.   The   equation   for 
mass  flow,  from  reference  11  is: 


w 


359    c    f    d2    Fa   Y    /hp  ,lb/hr 


where : 


c   Discharge  Coefficient 

f  =l/v/l-B!t     3=d/D 
d  and  D 


d 


T 


^F 


D 


^a   Thermal  Expansion  Factor 

y  (Ya  for  a  venturi)   Rdiabatic  Expansion  Factor 

h   The  pressure  irop  across  the  device,  in   inches   of 


water 


p   Density  of  the   fluid   being   metered,   pounds   p^r 
cubic  foot 


If  w1  is  the  entrained  flow  and  w2  the  primary,  thon  (w1 
uses  a  venturi,  w2  an  orifice) : 

r,  =  c,  f,  (dn)2  Fa-,  Ya   /h-, 
'2        52T2  ^2        F^2"^21  *2 


w- 
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Non-geometric  factors  for  the  primary  flow,  such  as 
discharge  coefficient,  could  be  reasonably  approximated  for 
all  tests  as: 


C2=.65 


Fa  2= 1.0 


Y2^.99 


In  the  case  of  the  entrained  flow  these  factors  had  to 
be  changed  for  each  test  axial  location,  due  to  the  widely 
varyinq  flow  requirements.  The  aeometric  factors  above  had 
the  followina  fixed  values: 


F1= 1.0  32  8   d1=.75 


F2=1.0792   d2=1.11 


Mow,    entrainment    ratio    is   commonly    written   as: 


=   w      +   w      =   W-, 


M 


0 


~w. 


v2 


Using  the  previous  formula  for   w1/w2, 
ratio  may  no*  be  written  as: 


M   =  A  /h 


M 


0 


IT 


1 


+1 


the   en tra inn  \ 


where   A   is   the   combination   of   the   ratios   of   factors 
previously   listed.   This  equation  was  used  to  calculate  ' 
nozzle  entrainment  ratios. 
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^PPEMHTX  C 


ENTRATNK2NT  CHAKBTP  DESIGN 


The  entrainment  chamber  design  was  closely  base-:!  on 
previous  work  by  Ricou  and  Spalding  (reference  7)  r  and 
Peschke  (reference  8).  Peschke's  report  especially  was  a 
good  source  of  specific  design  data,  while  Ricou  and 
Spalding  was  primarily  a  source  of  the  design  principles. 


The  en train me nt  chamber  consisted  essentially  of  a 
central  porous  cylinder,  closed  at  one  end  with  a  solid  wall 
through  which  the  nozzle  protruded,  the  other  end  being  open 
to  the  outside.  Orifica  plates  of  specific  diameters  were 
constructed  and  mated  to  cylinders  of  specific  length, 
length  and  diameter  specified  according  to  the  data 
available  in  Peschke.  These  orifice  cylinders  were  of  such 
outside  diameter  that  they  fit  closely  inside  the  porous 
cylinder,  into  which  they  were  inserted  for  the  measurement 
of  mass  flow  at  the  various  jet  axial  positions 
corresponding  to  the  orifice  diameters.  ^he  purpose  of  the 
orifice  was  to  seal  the  open  end  such  that  entrained  air 
could  only  enter  through  the  porous  cylinder.  The  purpose 
of  the  orifice  cylinder  was  to  seal  off  that  portion  of  the 
porous  cylinder  beyond  the  jet  axial  position  where  flow  was 
being  measured.  The  porous  cylinder  was  surrounded 
everywhere,  except  at  the  open  end,  by  a  box,  also  encasina 
the   nozzle  plenum,  which  served  as  the  membrane  air  supply. 

The  entrainment  chamber  was  about  4  by  U  feet  square  and 
H      feet  high,  supDorted  on  3  foot  legs.   The  porous  cylinder 
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was  10  inches  in  diameter  and  25  inches  long.  Orifice 
diameters  used,  corresponding  to  the  listed  jet  axial 
positions,  were  as  follows: 


Di  air.eter 


x/r> 


1.75 


2.5 


3.0 


5.0 


6.  o 


13.0 


The  entrainment  chamber  box,  legs,  orifice  elates  arid 
porous  cylinder  frame  work  were  all  constructed  of  wood. 
The  orifice  cylinders  were  made  of  fiberglass  sheet  and  the 
porous  covering  to  the  porous  cylinder  frame  work  was  tvo 
layers  of  dacron  crepe  cloth,  as  suggested  in  PeschkeTs 
report. 

One  critical  chamber  design  parameter  was  orifice 
diameter.  Tt  is  clear  that  too  small  an  orifice  would 
interfere  vi+h  the  jet  being  measured,  while  too  large  an 
orifice  would  not  properly  restrict  the  entrained  flow  to 
entering  only  through  the  membrane,  thus  preventing  the 
creation  of  a  readable  pressure  gradient  across  the  orifice. 
The  approach  to  this  problem  in  both  references  7  and  8  was 
to  find  an  optimal  orifice  size  a*  a  specific  jet  axijl 
location  by  testing  several  orifice  sizes.  To  minimize 
orifice  construction,  this  investigation  took  a  different 
approach.  That  was  to  use  the  reference  8  orifice  data  as  a 
starting  point,  and  vary  the  orifice  axial  location,  such 
that  a  readable  pressure  gradient  was  created  while  not 
interfering  with  the  jet.  In  both  approaches  the  attempt  is 
to  optimize  the  orifice  diameter,  jet.  axial  location 
combination . 
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One   primary   difference  between  this  chamber's  geometry 
ana  that  of  others  (references  7  and  8)  should  be  noted.   As 
shown   in   figure   8,   exit   orifices  were  positioned  in  ! 
porous  cylinder  at  various  distances  from  the  nozzle.     hat 
portion   of   the  porous  cylinder  extending  above  the  orifi 
plate  was  sealed  by  a  fiberglass  cylinder   attached   to   the 
orifice.    Thus   in   this  situation,  that  portion  of  the  jet 
beyond  the  orifice,  and  inside  the  porous  cylinder  exit,  •• 
not  in  a  truly  free  jet  flow  situation.   This  was  as  opoos^d 
to  the  reference  7   and   8   configurations   where   the   e 
orifice   corresponded   to   the   exit  of  ^he  porous  cylinder. 
Thus  in  those  experiments  the  jet  past  the  ori^ic?  was  in   a 
tru]y   free   jet   flow  situation.   The  configuration  usee?  in 
this   report   was   adopted   in   order   to    facilitate    the 
measurement   of   mass  flows  at  many  axial  locations.   Tt  was 
hoped  that  this  chamber's   geometry   would   nof   affect   The 
chamber   readings   greatly,   although   some  effect  was  to  be 
expected. 

An  important  proMes,  again  with  regard  to  the  chamber 
orifices,  concerned  the  measurement  of  the  pressure 
differencial  across  the  orifice.  Both  reference  7  and  8 
indicated  that  a  measurement  accuracy  of  .001  inches  of 
water  was  required.  Reference  7  used  a  micro  manometer 
while  reference  8  used  a  pressure  transducer  for  this 
purpose.  After  the  entrainmsnt  chamber  was  constructed 
experience  showed  that  such  accuracies  were  required.  This 
investigation  used  an  alchohol  micro  manometer  whi^h 
provided  an  accuracy  of  .000  1  inches  of  manometric  fluid. 
This  accuracy  was  obtained  through  the  combination  of  an 
inclined  tube,  hairline  indicator,  magnifying  lens  and  a 
micrometer  adjustment  to  indicator  level. 

Flow  measurement  devices  for  the  primary  and  entrain- 1 
flow  were  an  orifice  and  venturi  respectively.  These  were 
constructed  and  used  according   to   ASME   Power   Test   Codes 
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specifications     (reference     11).      Based     on    th( 
specifications,  an  inclined  water  manometer  was  pound  to   be 
necessarv   for   measuring  the  orifice  pressure  drop,  while  a 
vertical  mercury  raanemeter  was  sufficient  for  the  venturi. 
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APPENDIX  D 


PROGRAM  LISTINGS  AND  DOCUMENTATION 


This  appendix  contains  basic  program  and  input 
descriptions  for  the  two  programs  used,  from  references  1? 
and  13.  The  programs  also  are  listed  in  their  modi  Fled  form 
for  the  computation  of  a  ring  wing  at  angle  of  attach.  The 
proarara  statements  relevant  to  the  modification  are 
indicated  in  the  program  listings.  Also  indicated  in  the 
interior  flowfield  program  is  the  statement  uhose  omission 
caused  the  numerical  instability  noted  in  reference  12. 

The  inner  flow  field  program  consists  of  8  subroutines, 
and  a  main  program  which  manipulates  each  of  these.  The 
action  of  each  particular  subroutine  is  indicated  in  the 
program  listing.  Inputs  are  made  by  way  of  a  Fortran 
Namelist.   The  input  parameters  are,  in  proper  order: 

FSTRMN   Free  stream  roach  number. 

RFDFRQ    Ihe  reduced  frequency. 

RO    The  outer  cylinder  radius  to  length  ratio. 

RI    The  inner  cylinder  radius  to  length  ratio. 

NGRDFN    Grid  fineness,  or  the  number  of   grid   points 
taken  on  the  first  right  running  mach  line. 

N    Circumferential  mode  number. 

M    Axial  mode  number. 

H1    The  m  in  Qmn,  the  generalized  aerodynamic  force. 


m 


IPRTNT  An  output  parameter.  0  indicates  on]y  basic 
output.  1  indicates  output  including  such  things  as  Phi  an3 
its  derivatives  at  each  flow  field  point. 

Of  the  above  inputs,  in  the  modified  program,  the  inputs 

N,   M,   and   H1   must   be   set   to   1,   as   the   case  under 

consideration  is  not  panel  flutter  but  a  cylinder   at  angle 
of  attack.   The  other  inputs  are  arbitrary. 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  Q 

C  CYLINDRICAL  SHELL/PANEL  FLUTTER  PROGRAM 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  THIS  IS  THE  MAIN  PROGRAM 


C 


c 
c 


COMPLEX  PHK400)  ,DPHIS1(400)  fDPHIS2(400> 

DIMENSION  X<400) fR(400) fXPT(400> ,CPR(400) »CPI(400) 

DIMENSION  DUMMY (400) , DUNCE ( 4  00  ) 

CCMM0N/BLCK1/  FSTRMN t DELTAS , DX, DH, DSTSTR , RO  ,R 1 , REDFRO 

CCMMON/3LCK2/  Al , A2 , A3 , A4 , RNf  RM, RM1 »  FACT  I , F ACT2 

C0MM0N/BLCK3/  NGRDFN ,NPTS, LCCUNT , IHAVEP ,N,M,M1, IPRINT 

COMMON/ BLCK4/  PHI ,DPHIS1 »DPHIS2 

COMMON /BLCK5/  X, R , XPT  ,CPR, C PI 

DATA  PI/3.141593/ 

ENDPT(Y2tYl,XOfX2,DELX)  =  (  Y2  -  Yl  )*(  XO  -X2  ) /DELX 
+  Y2 


EPS    =    1  .E-04 
IER    =    0 

CALL     INPUT(IER) 
IF    (     IER     .EO.    1     )    GO    TO    12 
10  30    CONTINUE 
C 

IHAVEP    =    1 
LCCUNT    =    1 


C 

c 
c 


CALL  INITAL 
CONTINUE 
CALL  COMPXR 


C 

C      IS  THE  GRID  POINT  ON  THE  INITIAL  RIGHT  RUNNING  MACH 

C      LINE? 

C 

IF  ((  LCCUNT  .EQ.  1  )  .AND.  (IHAVEP  .LE.  NGRDFN  )} 
1    GO  TO  3 
C 
C      IS  THE  GRID  POINT  ON  THE  OUTER  CYLINDER  SURFACE? 


C 

C 
C 

c 

c 

c 

c 


IF  (  IHAVEP  .EQ.  0  I  GO  TO  5 
C 

IS  THE  GRID  POINT  ON  THE  INNER  CYLINDER  SURFACE? 

IF  (  IHAVEP  .EQ.  NGRDFN  )  GO  TO  6 
C 
C      THEN  THE  GRID  POINT  IS  IN  THE  GENERAL  FLOW  FIFLD. 


2  CALL  GENFPTC IER) 

IF  (  IER  .EO.  1  )  GO  TO  12 
IHAVEP  =  IHAVEP  ♦  1 
GO  TO  1 
C 

3  CALL  MACHLN 

IF  (  IHAVEP  .EO.  NGRDFN  )  GO  TO  4 

IHAVEP  =  IHAVEP  *■     1 

GO  TO  1 
C 
C      THIS  IS  THE  LAST  POINT  ON  THE  INITIAL  MACH  LINE. 


C 


LCCUNT  =  LCCUNT  +  1 
IHAVEP  =  0 
GC  TO  1 


5  CALL  RAD2 

C      IS  THIS  GRID  POINT  PAST  THE  FLEXIBLE  CYLINDER  LENGTH? 
C 
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c 
c 


c 
c 


c 
c 


IF  (  (X(  1) +  EPS)  .GE.  1.  )  GO  TO  8 
IHAVEP  =  IHAVEP  +  1 
GO  TO  1 

CALL  RAD1 

IHAVEP  =  0 

LCCUNT  =  LCOUNT  +  1 

GO  TO  1 

CONTINUE 

IF  (  ABSIXPT(LCOUNT)  -  1.  )  .LE.  EPS  )  GO  TO  10 

L  =  LCOUNT 

LM1  =  LCOUNT  -  1 

LM2  =  LCOUNT  -  2 

XFT(L)  =  1. 

CPR(L)  =  ENDPT(CPR(LM1)  ,CPR(LM2)  ,  1.  00  ,  XP  T<  LM1)  ,DSTSTP) 

CPI(L)  =  ENOPTCCPI (LM1) tCPI (LM2) » 1 . 00, X PT ( LM1 ) , DST STR < 


WRITE<6,105) 

DO  9  I  =  1,L 

ARG  =  RM1*PI*XPT( I ) 
C  C 

C      THE  FOLLOWING  TWO  CARDS  WERE  MODIFIED  FROM  THE        C 
C      ORIGINAL  REFERENCE  12  PROGRAM.  C 

C  C 

DUMMY (  I  )  =CPR ( I )* . 1*XPT( I ) 
9  DUNCE(I)=CPI( I)*.1*XPT(I) 


SUM1    =    0.5*<CPR(L)     *    CPP(LMl) )*( 1.00    -    XPT(LMl)} 
SUM2    =    0«5*(CPI(L)     *    CPKLM1)  )*<  1.00    -    XPT(LMl)) 

WRITE(6?U4)     FSTRMN, REDFRQ,RO,RI ,N 

CALL    QSF( DSTSTP  T  DUMMY , DUMMY , LM1 ) 
CALL    QSF(DSTSTR, DUNC E , DUNCE , LM1 I 
CREAL    =    0.5*(DUMMY(LM1)     *    SUM1) 
QIMAG    =    0.5=MDUNCE<LM1)    *    SUM2J 

W R I TE ( 6  v 1  I  5 )     M ,  M  i , QfcE  AL , M , M 1 , Q I M AG 
GO   TO    12 

10  CONTINUE 
XPT(L)     =    1. 
WRITE(6,105) 
DO    11     I    =     1, LCOUNT 
ARG    =    RM1*PI*XPT ( I ) 
DUMMY (  I  )=CPR(I  )*. l*XPT< I) 

11  DUNCE(  I)=CPI<  I)*.1*XPT(I1 

WRITE (6, 114)     FSTRMN, PEDF RQ, RO,RI ,N 

CALL    QSF( DSTSTR, DUMMY, DUMMY, L) 

CALL    QSF( DSTSTR, DUNCE, DUNCE ,  L) 

QREAL    =    0.5*DUMMY(LI 

QIMAG=     0.5*DUNCE(L) 

WPITE(6, 115)    M,M1,QREAL,M,M1,QIMAG 

12  CONTINUE 
WRITE(6,998) 

998  FORMAT( « 1» ,'CP    REAL') 

CALL     PLOTP(XPT,CPR, LCOUNT, 0) 
WRITE(6,999) 

999  FCRMATC ■ 1' ,'CP    IMAGINARY') 
CALL    PLOTP(XPT,CPI , LCOUNT, 0) 
STOP 

150    FCFMAT( F10.5)  .„„,.     ,,. 

105    FORMAT  (  •  1'  ,T19,M«  ,T27,'X»  ,T38,  «CPR«  ,T50  ,  '  CPI ■  ,//  ) 
110    FORMAT* «0' ,15X,I5,3(3X,F9.5) )  jt       ,    __ 

114    FORMAT(//»     ',T20,8MACH    NUMBER    =     ',F10.7,//'     ',120, 

1  'RPCUCED    FREQUENCY    =    «,F10.7,/'     ' ,T20 ,'( RADIUS/ • , 

2  'LENGTH)     OUTER    =    ',F10.7,/»     •,  T20  ,M  RACI  US/LENGTH  ) 
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3  ,»     INNER    =    ',F10.7,/'     '  , T20 , ' C  IRCUMF ERENTI AL    MODE', 

4  /»     • ,T34,« NUMBER    =     ',15,//) 

115    FOFMATCO'tTZOf'CllMZ.'.sn.MREAL    =     •  ,F  1  0.  5,  /  «  0  •  , 

1    T20, «0( » t  12 , •  ,«  ,12, • HMAG  =  ',F10.5,/) 
136  FORMAT(  •  1«) 
END 
C  THIS  IS  SUBROUTINE  INPUT 

C 

SLBROUTINE  INPUT(IER) 
C 

C      SUBROUTINE  INPUT  READS  INPUT  VARIABLES,  AND  DEFINES 
C      SEVERAL  CONSTANTS  USED  THROUGHOUT  THE  PROGRAM. 


C 


c 
c 

c 
c 


COMMON/ BLCK1/    FSTRMN, DELTAS ,DX, DH,DSTSTR ,RO,R I ,REDFRQ 
COMMON/ BLCK2/    Al , A2 , A3 , A4, RN,RM , RM 1 ,FACT lt  FACT2 
CCMMON/BLCK3/    NGRDFN , NPTS, LCCUNT , IHAVEP , N, M,M1< I  PRINT 
N  AMEL I  ST/ NAM  1 /    F  S  TRMN , REDFR  Q ,R  0 , R I , NGRDFN , N  ,  M , Ml , 
I  PR I  NT 

WRITE(6,80) 
READ(5,NAM1) 

WRITE(6,NAM1 ) 

FMANGL  =  ARSIN(l./FSTRMN) 
DH  =  (RO  -  RI)/FLOAT(NGRDFN) 

IF  (  DH  cLT.  0.  )  GO  TO  1 

DELTAS  =  DH*FSTRMN 
DX  =  DH/TAN( FMANGL) 
DSTSTR  =  2.*DX 
NFTS  =  1.  /DSTSTR  +  1 

IF  (  NPTS  .GT.  400  )  GO  TO  2 


BETA  =  SQRT( FSTRMN-FSTRMN  -  1.} 

RM  =  FLCAT(N) 

RN  =  FLOAT (M) 

Rf'l  =  FL0AT(M1  ) 

FACTi  =  REDFRQ*REDFRQ 

FACT2  =  RN*RN/(FSTRMN*FSTRMNI 
C 

C      COMPUTE  CONSTANTS  TO  BE  USED  IN  THE  COMPUTATIONAL 
C      MOLECULES. 
C 

Al  =  0.2  5*DELTAS/FSTRMN 

A2  =  0.5*FSTRMN/BETA 

A3  =  REDFRG*DELTAS*A2 

A4  =  0.5-FSTRMN 

RETURN 
C 

1  IFR  =  1 
WRITE(6,100) 
GO  TO  3 

C 

2  IER  =  1 
WPITE(6,110) 

3  RETURN 

100  FORMAT<//»  •  ,35X ,  « ABNORMAL  TERMINATION:  INPUT  RADII  ', 

1101F0RMAT( //■  •  ,35X, •ABNORMAL  TERMINATION:  MUST  INCREAS1, 
1    »E  THE  SIZE  OF  VECTORS  XPT,  CPR,  AND  CPIS//> 

END 
C  THIS  IS  SUBROUTINE  INITAL 

SLBROUTINE  INITAL 

C      SLBROUTINE  INITAL  INITIALIZES  ALL  FLOW  FIELC 
C      QUANTITIES. 

C 

CCMPLEX    DPHIDR,DPHIDX 
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CCMPL5X    PHK^.0  0)  , DPHI S 1 ( 400 ) , D PH  I  S2 ( 400 ) 

DIMENSION    X(400) ,R<400) ,XPT<400) ,CPR(400 ),CPI(400) 

S£'SH2tK  £.L£K  l  '    F  S  TRMN ,  DE  LT  AS  ,  DX  ,  DH  ,  DSTS  T  R  ,  R  0  ,  R  I ,  k  E  DF  r 

CCMM0N/BLCK2/    AI , A2 , A3 , A4, RN,RM » RM1 , FACT  1 ,FACT2 

£S  fJ2^/,BLCK3/    ^C-RDFN,NPTS,LCQUNT,  IHAVEP , N, M  , Ml ,  IPRINT 

CCMM0N/BLCK4/     PHI ,DPHI SI , DPHIS2 

COMMON/ B LC K5/    X , R , XP  T  ,  CPR , CP I 

DATA    PI/3.141593/ 
C 

X( 1}    =    0. 

R(l )    =    RO 

XFT(l)     =    0. 
C 

C      ALONG  THE  INITIAL  RIGHT  RUNNING  CHARACTERISTIC  PHI  AND 
C      DPHIS2  ARE  ZERO. 
C 

PHH1)  =  (0.,0.) 

DPhIS2( 1  )  =  ( 0.,0. ) 
C 
C 

AA=.2/FSTRMN 

DPHISK  1)     =    CMPLXt  AA,0.  ) 

DPHIOk  =  A4*DPHIS1<1) 

DPHIDX  =  A2*DPHIS1U ) 

C  PR  i  1  )  =  -2.*  RE  AL ( D  PHI DX ) 

CPKU  =  -2.*AIMAG(DPHIDX> 
C 

IF  {  IPRINT  ,NE.  1  )  RETURN 

WRITE (6, 90  J 

WRITE (6, 95) 

WRITE  (6,  100)  X(  1.)  v  DPHISK  i)  ?DPHIS2(  I)  ,  PHI  (  1), 
1    R( 1) , DPHIDR,  DPHIDX 

WRITE<6,110)  CPR(1),CPI(1) 

RETURN 

90  FGBMATC  I1  ,T10,'X«  tT22»,0(PHI  )/D<Sl  )  «  fT46t 

1    «0( PHI)/D(S2J  «  ,T70,«PHI!  ) 
95  FORMAT ( •  «  , T 10 , * R •  ,T23 , • D ( PHI ) /DR« ?T47t ■  D< PHI  I /OX1  ,  /  ) 
100  FORMATOH  I  , 2X , F6 .4t4X,3 ( F 12. 5 »F9.5  ,3X )  i  / '  ,,4X,F8,4, 

1    4Xt3(F12.5,F9.5,3X) ) 
110  FORMAT ( »0' ,10X»«  INITIAL  VALUE:   INNER  SURFACE  OF • t 

1  ■  OUTER  CYLINDERS/'  '  t30X  ,  •  CP(  RE  AL )  =  »,F8.4,10X, 

2  CCP{  IMAG)  =  •  ,F8.4,//) 
END 

C  THIS  IS  SUBROUTINE  COMPXR 

C 

SUBROUTINE  COMPXR 
C 

C      SUBROUTINE  COMPXR  COMPUTES  THE  (X,R)  COORDINATES  OF  A 
C      GRID  POINT. 
C 

DIMENSION  X(400  )  tR(400) ,XPT(400  )  ,CPR< 400 ) f CPI ( 400) 

CCMM0N/BLCK1/  FSTRMN, DELTAS, OX, DHfDSTSTR , RO,RI t  REDFRC 

CCM.M0N/BLCK3/  NGRDFN  ,NPTS,  LCOUNT  ,  IHAVEP  ,  N,  M  ,  Ml ,  I  PR  INT 

CCMMQN/BLCK5/  X, R , XPT  ,CPR , CP  I 
C 

I  =  IHAVEP  *  1 

J  *  IHAVEP 
C 

C      IF  THIS  IS  THE  FIRST  POINT  ON  A  NEW  RIGHT  RUNNING  MACH 
C      LINE,  GO  TO  1. 


C 


IF  (  J  .EQ.  0  )  GO  TO  1 

X( I)  =  X( J)  +  DX 
R( I)  =  R( J)  -OH 

RETURN 

R(I)  ~  RO 

X( 1)  =  X( 1)  +  DSTSTR 

RETURN 

END 


THIS  IS  SUBROUTINE  MACHLN 
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c 

c 
c 
c 
c 


c: 

c 


c 
c 
c 

c 

c 
c 
c 
c 
c 
c 
c 


c 
c 

c 

c 


SUBROUTINE  MACHLN 

SUBROUTINE  MACHLN  COMPUTES  THE  VALUES  0^ 
FIELD  QUANTITIES  ALONG  THE  INITIAL  RIGHT- 
LINE  WHERE  PHI  AND  DPHIS2  ARE  ZERO. 


THE  FLOW 
•RUNNING  MACH 


CCMPL 
COMPL 

CCMPL 
D I  ME  N 
CCMMC 

cor 

CCMMG 
CCMMO 
COMMO 


EX  A,C 
EX  DPH 
EX  PHI 
SIGN  X 
M/BLCK 
N/BLCK 
N/BLCK 
N/BLCK 
N/BLCK 


IDR, DPHIDX 

(400)  , DPHIS1 (400) ,D PHIS 2(400) 

(400)  ,R(400),XPT(400)  «CPRU00  KCPK400  ) 

1 /  FS T RMNt OE L TAS , DX, DH, DSTST  R , RO ,RI , REDF  RC 

2/  Ai , A?, A3, A4,RN,RM,RM1 ,FACT1, FACT2 

3/  NGRDFN,NPTS,LCGUNT  ,  IHAVEP , N, M , Ml ,  I  PRINT 

4/  PHI ,DPHIS1 ,DPHIS2 

5/  X,R,XPT,CPR,CPI 


I  =  IHAVEP  4-  1 
J  =  IHAVEP 

CALCULATE  THE  COEFFICIENTS  ASSOCIATED  WITH  THE  RIGHT- 
RUNNING  MACH  LINE. 
LINE. 

AA1  =  1.  -  Al/RU) 

A  =  C 4PLX( AA1, A3) 

CI  =  1.  *  Al/RU  ) 

C  =  DPHIS1 (J)*CMPLX(Clt-A3) 

Phim  =  (0.,0.) 
DPHIS2(I)  =  (0.,0.) 
DPHISK  I  )  =  C/A 

DFHIDR  =   A4*DPHIS1(I ) 
DPHIDX  =  A2*DPHIS1(I  ) 

IF  (  I  PRINT  .NE.  1  )  RETURN 

WRITE  (6,  100)  X(  I)  ,  DPHISK  I  ),  DP  HIS2(  I), PHH  I), 
1    R( I ) , DPHIDR, DPHIDX 
WRITE(6,111J 
RETURN 
130  FCRMAK3H  M  t  2X,  F  8  .4,  4X,3  (  F  12.  5,  F9  .  5  ,3X  )  ,  /  »  »,4X,F8.4, 

1    4X,3(F12.5»F9.5,3X) ) 
111  FCPMATC '0*  ) 
END 

THIS  IS  SUBROUTINE  RAD1 

THIS  IS  SUBROUTINE  RAD1 
SUBROUTINE  RAD1 

SUBROUTINE  P.AD1  COMPUTES  THE  FLOW  FIELD  QUANTITIES  AT 
THE  SURFACE  OF  THE  INNER  RADIUS  WHERE  THE  BCUNDARY 
CONDITION  IS   DEPENDENT  ON  N  ,  THE  CIRCUMFERENTIAL 
MODE  NUMBER.   IF   N   IS: 

EVEN   —   DPHIDR  =  0,  AND  THUS  DPHIS1  =  DPHIS2. 

ODD    —   CP  =  0. 

CCMPL EX  A,B.C,TANR,DPHIDR,DPHI DX , DELTA, G AMMA, IK 
CCMPL  EX  PHI  (400)  ,  DPHISK  400)  ,  DPH  IS  2(4001 
DIMENSl ON  X(400) , R (400 1 , XPT ( 400 ) ,CPR(400l ,CPI(400) 
CCKM0N/BLCK1/  FST RMN, DELTAS , DX, DH, DSTSTR , RO, RI , REDFRC 

A  i , A  2 , A3 , A4 , RN, P  M , RM 1 , F ACT  1 r  F  ACT2 
NGRDFN,NPTS,LCCUNT, IHAV EP » N, M , Ml , I  PRINT 
PHI ,DPHIS1 ,DPHIS2 
X,R,XPT,CPR,CPI 


CCMM0N/BLCK2/ 

CCMM0N/BLCK2/ 
CCMM3N/8LCK4/ 
C0MM0N/BLCK5/ 


I  =  IHAVEP  +  1 
J  =  IHAVEP 

CALCULATE  THE  COEFFICIENTS  ASSOCIATED  WITH  THE  RIGHT- 
RUNNING  MACH  LINE. 
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c 
c 
c 
c 
c 
c 


c 
c 
c 
c 
c 
c 


c 
c 
c 


c 
c 
c 


c 

c 

c 


c 
c 

c 
c 
c 
c 
c 


THE    FOLLOWING    VARIABLE,     IK,     WAS    FOUND    TO    BE 
UNDEFINED    IN    THE    ORIGINAL    REFERENCE    12     PROGRAM,     A 
WAS    THE     REASON    FOR    A    NUMERICAL     INSTABILITY    AFTER 
FIRST    REFLECTION. 


ND 
THE 


IK=CMPLX(0.0,REDFRQ) 

AA1    =     1  .    -    Al/R< I ) 

A    =    CMPLXC AAl, A3) 

Bl    =    A1/RU)    -    0.25*DELTAS*DELTAS*(FACT1    - 
1  FACT2/(R(  I)*R(  I))  ) 

B    =    CMPLX(B1,A3) 

CI    =    1.     <r    A1/R(J) 

C2    =    -A1/R{J)     <-    0.25*DELTAS*DELTAS*(FACT1    - 
1  FACT2/( R(  I)*R<  I  )  ) ) 

C3    =    (2.*FACT1    -    FACT2/(R( J)*R( J) )    - 
1  FACT2/(R(  I)*R(  I  )  )  )*DELTAS*0.5 

C  =  OPHIS1 ( J)*CMPLX(C1,-A3)  +  DPHIS2 { J  I *CMPLX( C2 , -  A3 ) 
1    +  PHIl J)*CMPLX(C3,0.) 


IF  (  MOD(N,2J  ,EQ.  0  )  GO  TO  1 

APPLY  THE  BOUNDARY  CONDITION  FOR  ODD  VALUES  OF  N. 
CP  =  (IK) PHI  +  DPHIDX  =  0.   USE  THIS  B.C.,  THE  FI 
DIFFERENCE  EQN.  FOR  DPHIDX,  AND  THE  EQN.  FOR  THE 
RUNNING  MACHLINE  TO  SOLVE  FOR  DPHIS1  AND  DPHIS2. 

DELTA  =  A2  +  0 . 5*1 K*DELTA S 

GAMMA  =  IK*(  PHI (J)  *  0.5*DELTAS*DPHIS?( J)  ) 


NITE 

RIGHT 


SOLVING  FOR  DPMI  SI  AND  DPHIS2 


DPHISK  I  )  =  ( 
DPHIS2( II  =  ( 


C  +  B*GAMMA/l)ELTA  )/( 
C  -  A*DPHIS1(I)  )/B 


A   -     B*A2/DELTA     ) 


DPHIDX  =  A2*(  DPHISKI)  *■  DPHIS2(I)  ) 
DPHIDR  =  A4*(  DPHISKI)  -  DPHIS2(I)  ) 
GC    TO    2 

APPLY    THE    BOUNDARY    CONDITION    FOR    EVEN    VALUES    OF 


N. 


1  D PHIS  It  I )  =  C/(A+B) 
DPHIS2(  I)  =  DPHIS1 (I  ) 
DPHIDR  =  (0..0.) 

CPHIDX  =  A  2*  (DPHISK  I  )  *■    DPHIS2(I)) 

PHI. 

2  PHKI)  =  PHI(J)  <-  0.5*(DPHIS2(  J)  +  DPH  I  S  2(  I  )  )  *DELT  AS 

IF  (  IPRINT  .NE.  1  I  RETURN 

WRITE  (6, 100)  X(  I  )  ,  DPHISK  I  )  ,DPHIS2(I)  ,  PH  I  (  I  ), 
1    R( I )  , DPHIDR, DPHIDX 
RETURN 
100  F0RMAT(3H  RI , 2X, F8 .4, 4X , 3 ( F 1 2 . 5 , F9 . 5 , 3X ),/ '  ',4X,F8.4, 
1    4X,3(F12.5,F9.5,3X) ) 
END 

THIS  IS  SUBROUTINE  RAD2 

SUBROUTINE  RAD2 

SUBROUTINE  RAD2  COMPUTES  THE  FLOW  FIELD  QUANTITIES  AT 
THE  INNER  SURFACE  OF  THE  OUTER  CYLINDER,  WHERE  THE 
FLOW  TANGENCY  CONDITION  PRESCRIBES  DPHIDR. 

COMPLEX  D,E,F,TAMR 

CCMPLEX  DPHIDR, DPHIDX 

COMPLEX  PHH400)  ,  DPHISK  400)  ,DPHIS2(  400) 

DIMENSION  X(400) , R ( 400) , XPT ( 400 ) ,CPR( 400 ) ,C PI ( 400 ) 

CCMMPN/3LCK1/  FSTRMN, DELTAS, DX , DH,DSTSTR ,R0 ,R I , RE DFRQ 

CCMUCN/BLCK2/  Al , A2 , A3 , A4 , RN, RM , RMl , F ACT  1 , F ACT2 

COMMON/ BLCK3/  NGRDFN , NPTS , LCOUNT  ,  IHAVEP , N, M  ,M1 , 1 PRI NT 
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CCMM0N/BLCK4/  PH I , DPH I  SI , DPH  I S2 

C0MM0N/BLCK5/  X, R , XPT , CPR » C P I 

DATA  PI/3.141593/ 
C 

I  =  IHAVEP  4-  1 

J  =  IHAVEP  +■  2 

K  =  LCOUNT 

ARG=RM*PI#X< I) 
C 
C 

C      CALCULATE  THE  COEFFICIENTS  ASSOCIATED  WITH  THE  LEFT 
C      RUNNING  MACH  LINE. 
C 

Dl  =  A1/R(I)  *  0.25*DELTAS*DELTAS*(FACTl  -  FACT2/(R(I) 
I    *  R  (  I  }  )  ) 

D  =  CMPLX(01,-A3J 

El  =  -1.  -  Al/R( I ) 

E  =  CMPLX(E1 t-A3) 

Fl  =  -Al/RU)  ~  0.25*DELTAS*DELTAS*(FACT1  -  FACT?./ 
1    (R(I)*RU))) 

F2  =  -1  .  4-  Al/R(  J  ) 

F3  =  (2.*FACT1  -  FACT2/(R(  I  )  *R  (  I  )  )  -  FACT2/(R(J)« 
I    R(J) ))*DELTAS*0.5 

F  =  DPHIS1(J)*CMPLX(F1,A3)  ♦  DPH I S2{ J) *C MPLX( F2 , A3 i     - 
1    PHI ( JJ*CMPLX(F3iO. J 
C 

C      CALCULATE  THE  FLOW  TANGENCY  CONDITION   -   T ANR 
C 

C      BOUNDARY  CONDITIONS  FOR  PANEL  FLUTTER 
C 

C  C 

C      THE  FOLLOWING  THREE  CARDS  WERE  MODIFIED  FROM  THE      C 
C      ORIGINAL  REFERENCE  12  PROGRAM.  C 

C  C 

ALPHAS  1 

TAKR1=ALPHA 

TANR2=REDFRQ*X(I )* ALPHA 

TANR  =  CMPLX(TANR1»TANR2) 
C 

DPH  I  SKI)    =    (    F    4-    TANR*E/A4    )/(     D    +    E    ) 

DPH1S2U)    =    (     F    -    TANR*D/A4    )/<     D    4-    E    ) 

DPHIDR    =    A4*(DPHIS1 ( I )    -    DPHIS2(I>) 

DPHIDX    =    A2*(DPHIS1(I)     *    DPHI.S2U)) 
C 

C      INTEGRATE  ALONG  THE  LEFT  RUNNING  MACH  LINE  TO  FIND  PHI 
C 

PHI(I)  =  PHI(J)  4-  0.5*(DPHIS1(J)  4  DPHISM  I  ))*DELTAS 
C 
C      CALCULATE  THE  PRESSURE  COEFFICIENT  AT  THIS  GRID  POINT. 


C 


c 


yrT|  1/  I   =   y  /  T  j 

CPR(K)     =    -2.*<-REDFRQ*AIMAG(PHI( I M    +    RE AL ( CPHIDX) ) 
CPI(K)     =    -2.*(REDFRQ#REAL(PHI(I) )     4-    AIMAGt DPHIDX) ) 


IF    (     IPRINT    .ME.     1     )     RETURN 

WRITE  (6,100)     X(  I)  ,  DPH  I  SKI  )  ,DPHIS2(I)  ,  PHKI  ), 
1  R( I ) ,DPHIDR, DPHIDX, TANR 

WRITE(6,110)    CPR(K)tCPKK) 
RETURN 
100    FORMATOH    RO,  2X  ,  F8  .4,  4X,  3  <  Fl  2  .  5  ,  F9.5  ,  3X  )  ,  /  ■     *,4X,F8.4, 

1  4X,3(F12.5,F9.5,3X) ) 

110    FCRMAT(  «0»  ,10X,» INNER    SURFACE    OF    CUTER    CYLINDER', /'     »  , 
1  30X,  'CP(REAL)     =     '  ,F3.4,10X, »CP(  IMAGI     =    ',F8.4,//) 

END 
C  THIS    IS    SUBROUTINE    GENFPT 

C 

SLBPOUTINE  GENFPT( IERJ 

C      SLBROUTINE  GENFPT  COMPUTES  THE  FLOW  FIELD  QUANTITIES 
C      AT  A  GENERAL  FLOW  FIELD  POINT. 


CCf'PLEX  PHIL 
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c 

c 


c 
c 
c 
c 


c 
c 
c 
c 
c 


c 

c 
c 
c 


CCFPLEX    A?B,C,D,E,F, DENOM 

CCMPLEX    DPHIDP .DPHIDX 

COMPLEX    PHH400)  ,DPHI  SU400)  ,DPHIS2<400) 

DIMENSION  X(400l , R<400) , XPT { 40  0 ) ,CPR I 400 ) , C PI ( 400) 

CCMMON/BLCKl/  FS  T  RMN  t  DELTAS  ,  OX  ,  DH,  OSTST  R  ,  RO  ,  R  I  ,  REDFRC 

Al,A2tA3,A4fRN,RM,RMl,FACTl,FACT2 

NGRDFN,NPTS,L COUNT  ,  IHAVEP  ,  N,  M  f  M  J. ,  I  PR  INT 

PHI ,DPHISl ,DPHIS2 

X,R,XPT,CPR,CPI 


COMMON/ BLCK 2/ 
CCMM0N/BLCK3/ 

CCMMCN/BLCK4/ 
CCM-10N/BLCK5/ 


I    =    IHAVEP    * 
J    -    IHAVEP 
K    =    IHAVEP    + 


1 


COMPUTE    THE    COEFFICIENTS    OF    OPHISI    AND    DPHIS2. 
F£CT3    =     FACT 2/ { R(  I  )*R( I)  ) 
FACT4    =     FACT2/(R(K)*R(K) ) 
FACT5    =    FACT2/(R< J)*R< J) ) 

ATB    AND   C    ARE    ASSOCIATED    WITH    THE    RIGHT     RUNNING    MACH 
LINE,    WHILE    D,E    AND    F    ARE    ASSOCIATED    WITH    THE    LEFT. 

AA1    =    1.  .    -    Al/R(  I  ) 
A    =    CMPLX1 AAlt A3) 

Bl    =    -AA1    *    1.     -    0.25*0ELTAS*DELTAS*(FACT1     -    FACT3) 
B    =    CMPLX(B1,A3) 
Ci    =    1.     4-    Al/P.U) 

C2    =    -CI    *    I  .    <-    0.25*DELTAS*DELTAS*(FACT1    -    FACT3) 
C3    =    (2.*FACT1    -    FACT3   -    FACT5 )*DELTAS*0 .5 
C    =    DPHISll'J)*CMPLX(Ci,-A3)     *    DPHIS2(  J)  *CMPLX(C2,-A3) 
1  +    PHI(  JKCMPLX(C3,0.  ) 

Dl    =    Al/RCI)    +-    0.25*DELTAS*DELTAS*(FACT1    -     FACT3) 
D    =    CMPLX(D1,-A3 ) 
El    =    -1 .    -    Al/R( I) 
E    =    CMPLX! El, -A3  * 

Fi    =    -A1/R(K)    -    0.25*DELTAS*DELTAS*(FACT1    -    FACT3) 
F2    =    -1.     <■    Al/MK) 

F3    =    (2.*FACTl    -    FACT3    -    FACT4)*DELTAS*0 .5 
F    =    DPHISi(K)*CMPLX(Fl,A3)     +    D PHI S2( K)*CMPLX( F2, A3)    - 
1  PHI(K)*CMPLX(F3,0. ) 

USE    CRAMERS    RULE    TO    SOLVE    FOR     DPHISi    AND    DPHIS2 
PROVIDED    THAT    THE    DETERMINANT    OF    COEFFICIENTS    IS    NCN- 
SINGULAR. 


DEKCM  =  A*E  -  D*B 

IF  (  CABS(  DENOM  )  .LT.  l.E-05  )  GO 
DPHISI (IJ  =  (C*E  -  B*F) /DENOM 
DPHIS2CI)  =  <A*F  -  C*D  ) /DENOM 
DFHIDR  =  A4*(DPHISim  -  DPHIS2U)) 
CPHIDX  =  A2*(DPHIS1<  I  )  +•  DPHIS2(I)) 


TO    I 


INTEGRATE    ALONG 
PHI. 


THE    RIGHT    RUNNING    MACH    LINE    TO    FIND 


PHKI)     =    PHI(J)     *•    0.5*(DPHIS2( I )     +    DPHI  S  2(  J  )  )  *DE  I.TAS 
PHIL    =    PHI(K)     «■    0,5*(DPHIS1(K)     +    DPHISI  (  I  ))  *DELTAS 

IF    (     IPRINT     .ME.     1     )     RETURN 

WRITE ( 6, i 00)     X( I ) ,DPHIS1( I ) ,DPHIS2( I),PHI( I ), 
i  R(  I)  ,DPHIDR,DPHIDX,  PHIL 

WPITE<6,111) 
RETURN 


100 
110 


1  IER  =  1 

WRITE(6,110)  DENOM 
RETURN 

F0RMATC3H  G  ,  2X , F 8  .4, AX , 3 ( F 1 2. 5 , F9 . 5 , 3X  )  , 
1    4X,3(F1 2. 5,F9.5,3X) ) 


/'  »  T4X,F8.4, 


FORMAT ( «0«  ,///,«  '  ,35X,1 ABNORMAL  TERMINATION 


•DETERMINANT 


OF'COEFFICIENTS 


IS  SINGULARS/ 


i 

0« 


50  X, 
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2    'DENOMINATOR  =  • ,E12.5) 
111  FORMAT ( •O' ) 
END 
C  THIS  IS  SUBROUTINE  CSF 

C         SUBROUTINE  QSF 

SUBROUTINE  0 SF ( H , Y , Z , ND I M ) 
C  CALL  QSF  <H,Y,Z,NDIM) 

C         DESCRIPTION  OF  PARAMETERS 

C  H         THE  INCREMENT  OF  ARGUMENT  VALUES. 

C  Y         THE  INPUT  VECTOR  OF  FUNCTION  VALUES. 

C  Z         THE  RESULTING  VECTOR  OF  INTEGRAL 

C  VALUES.   Z   MAY  BE  IDENTICAL  WITH   Y. 

C  IDENTICAL  WITH  Y. 

C  NDIM    -  THE  DIMENSION  OF  VECTORS  Y  AND  Z. 

C 

C         REMARKS 

C  NO  ACTION  IN  CASE  NDIM  LESS  THAN  3. 

C 

c 
c 

DIMENSION  Y( 1) ,Z(1 ) 
C 

HT=.3333333*H 

IF(NDIM-5) 7T8, 1 
C 

C      NDIM  IS  GREATER  THAN  5.  PREPARATIONS  OF  INTEGRATION 
C      LCCP 
C 

1  SUM1=Y(2)*-Y(2) 
SUMlsSUMl+SUMl 
SUM1=HT*(Y(1)4-SUM1«-Y(3)  ) 
AUX1=Y(4)*-Y(4) 
AUX1=AUX1+AUX1 
AUX1=SUN;JL*HT*<Y<3)+AUX1+Y<5)  J 

AUX2=HT*(Y<  1  )<-3.  87  5*1  Y(2)+Y(5)  J  +2 .625*( Y < 3 } +Y< 4)  )  «■ 
1    Y(&)) 

SUf2=Y<5»+Y(5) 

SUM2=SUM2+SUM2 

SUM2=AUX2-HT*(Y< 4) +SUM2+Y( 6  ) ) 

Z(1)=0- 

AUX=Y(3)+Y<3) 

AUX=AUX+AUX 

Z  (  2 )  =SUM2-HT*  (  Y (  2  )  4-AUX+Y ( 4  )  ) 

Z(3)=SUM1 

Z(4)=SUM2 

IF(NDIM-6)5,5,2 
C 
C      INTEGRATION  LOOP 

2  DO  4  I  =  7,  NDIM,  2 
SUM1=AUX1 
SUM2=AUX2 

AUX1-YC I-1H-YC 1-1) 
AUX1=AUXH-AUX1 

AUXi  =  SUMH-HT~-(Y(  I-2)«-AUXH-Y(I  )  ) 

Z(  1-2) = SUM  I 

IF( I-NDIM)3,6,6 

3  AUX2=Y( I )+Y( I) 
AUX2=AUX2+AUX2 
AUX2=SUM2«-HT*(Y(  1-1)4- AUX2+Y<  1*1)  ) 

4  Z(  l-l)=SUM2 

5  ZCNDIM-1 )=AUX1 
Z(NDIM)=AUX2 
RETURN 

6  Z(NDIM-1)=SUM2 
Z(NDIM)=AUX1 
RETURN 

C      END  OF  INTEGRATION  LOOP 
C 

7  IFCNDIM-3)  12, 1 1 1 8 
C 

C  NDIM    IS    EQUAL    TO    4    OR    5  , 

8  SUM2=1.125*HT*(Y(1 ) +Y( 2 ) +Y ( 2 } +  Y<  2 ) +Y( 3 ) 4- Y( 3  )+Y( 3)  4- 
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1  Y(4J) 

SUM1=Y(2)+Y<2) 

SUfn=SUMl  +  SUMl 
SUM1  =  HT*<  YU)+SUNIH-Y< 3) ) 
Z(1)=0. 

AUX1=Y( 3) *Y(3) 
AUX1=AUXH-AUX1 

Z<2)=SUM2-HT*(Y<  2  )  +  AUX  1*Y(  4  )  ) 
IF(NDIM-5)  10,9,9 
9    ALX1  =  Y(A)«-Y(4) 
AUX1=AUX1+AUX1 
Z { 5 l  =  SU M i+HT*  <  Y  <  3 ) +AUX 1+Y  <  5  ) ) 

10  Z(3)=SU"1 
Z(4)=SUI 
RETURN 

C 

C  NDIM    IS    EQUAL    TO    3 

11  SUM1=HT*( 1.25*Y( l)tY(2)tY( 2$-.25*Y(3) ) 
SUM2=Y( 2)+Y(2) 

SUM2=SUM2+SUM2 
Z<3»=HT*(Y(1  KSUM2+Y(3)  ) 
Z(  1)=0. 
Z(2)=SUM1 

12  RETURN 
END 
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The  outer  flowfield  program,  from  reference  13 ,  consists 
of  5  subroutines  and  a  main  program  to  manipulate  the 
subroutines.  The  action  of  each  oarticular  subroutine  is 
indicated  in  the  program  listing.  Inputs  reguircd  to  the 
program  must  be  on  two  cards,  and  are  as  follov.'s: 

DATE 

HACH   Mach  number. 

RADIUS    Radius  to  length  ratio. 

K    Reduced  frequency. 

MFREQ    Axial  mode  number. 

UP.     The  r  in  Qmr,  the  generalized  aerodynamic  force. 

NFREQ    Circumferential  mode  number. 

FINGRD    Grid  fineness,  essentially  the  number  of  err  id 
points  along  the  cylinder  surface. 

The  run  date  is  typed  onto  the  first  card.  The  second  card 
contains  the  above  parameters  in  the  following  format: 
F20.8,  2^10.3,  415.  Of  the  above  parameters,  in  the 
modified  program,  MFREQ,  MR,  end  NFREQ  must  be  set  to  1,  as 
the  case  under  consideration  is  not  panel  flutter  but  a 
cylinder  at  angle  of  attack.  Other  parameters  are 
arbitrary. 
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